NRQCD and Static Systems - A General Variational Approach * 
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We present initial results from Monte Carlo simulations of NRQCD-light, static-light, and NRQCD-NRQCD 
mesons, using a variational technique (MOST), as part of our ongoing calculation of the /b decay constant. 
The basis states for the variational calculation are quark-antiquark operators separated by all possible relative 
distances not equivalent under the cubic group (for example, for a 20^ lattice there are 286 operators). The efficacy 
of the method is demonstrated by the good plateaus obtained for the ground state and the clean extraction of 
the wave functions of the ground and first radially excited state. 



Theoretical uncertainties are starting to dom- 
inate the determination from experiment of the 
Cabibbo-Kobayashi-Maskawa quark mixing ma- 
trix, with a major bottleneck being the reliable 
calculation of the heavy-Hght pseudoscalar de- 
cay constant, /s [1]. An essential part of its 
determination on the lattice is its value in the 
static [2] and nonrelativistic (NRQCD) [3] ap- 
proximations. Unfortunately, the static approxi- 
mation is plagued by the particularly bad signal- 
to-noise ratio for the correlation functions. The 
use of several standard lattice techniques has led 
to differing results from various groups in report- 
ing values of the heavy-light pseudoscalar decay 
constant in the static approximation on the lat- 
tice [4]. This can be traced to the difficulty in 
isolating unambiguously the ground state at early 
Euclidean time separations [5]. 

The FNAL group vi^ere able to get reliable re- 
sults for fg^^ [6] from a variational calculation 
using a basis set of trial smearing functions de- 
signed to have good overlap with the physical 
states, having been constructed as wave functions 
of a simple relativistic quark model Hamiltonian. 

The present authors used a complementary and 
more general approach [7]: "MOST" (Maximal 
Operator Smearing Technique). With MOST, op- 
timally smeared operators are generated without 
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any ansatz, and emerge naturally and directly 
from the QCD dynamics. This variational cal- 
culation uses a complete basis set of Qq interpo- 
lating operators 



*Si (x, t) 



Y,Si{r)Q{x + r,t)Tq{x,t) (1) 



Each Si is a sum of delta functions of relative 
separations r equivalent under cubic symmetry. 
The complete set of these "onion shells" is used 
to construct the N x N two-point correlation ma- 
trix C(i) (with TV = 165 and 286 for lattice vol- 
umes of 16^ and 20^ respectively). This matrix 
contains all of the information available from two- 
quark interpolating fields. It can be analyzed at 
leisure on a workstation using a variety of dif- 
ferent techniques including, for example, multi- 
exponential fits to a submatrix defined by model 
wave functions, or an iterative method based on 
using a measured wave function as a source. The 
best method will be a compromise between max- 
imizing the information obtained from the ma- 
trix and minimizing the effects of the noisy el- 
ements. We use a general approach to extract 
the amplitudes Am = (vacl^^Jn) (which are the 
desired Bethe-Salpeter wave functions). If the 
wave functions were orthogonal they could be ex- 
tracted from a straightforward diagonalization of 
the matrix; however, since they are not, the dual 
wave functions, Z, are required. These satisfy 
Z^A = 1, and can be obtained by diagonaliza- 
tion of C(io)"^/^C(i)C(io)"^/^[8,9]. These "S"- 
sources - optimally-smeared interpolating oper- 



ators, are then used to construct local-smeared 
and smeared-smeared correlators, from which Jb 
is extracted in the usual way. The optimal inter- 
polating operators project out ground and excited 
states at early Euclidean times when statistical 
noise is smallest. 

The method can be appHed to create optimal 
smearing functions for lattice simulations of any 
heavy-quark effective theory, whenever it is com- 
putationally inexpensive to compute many heavy 
quark propagators. It is then feasible to use 
the complete basis, using fast-Fourier transforms 
to evaluate the convolutions. Here we focus on 
the application of MOST to heavy-Hght mesons, 
where the heavy quark is described by a nonrela- 
tivistic action (NRQCD) [3,10], and whose prop- 
agation can be computed at low cost, relative to 
that of a light Wilson quark, by recursion [3]. 
As a by-product, one can construct NRQCD- 
NRQCD mesons for free. 

Light quark propagators were computed for 
K = 0.154, 0.155, and 0.156 with the Wilson 
fermionic action in a background of 32 quenched 
(Coulomb gauge-fixed) configurations (selected 
every 1000 pseudo-heat-bath sweeps after 5000 
thermalization sweeps) with (3 = 6.0 on a 20^ x 30 
lattice. 

Reference [7] contains our preliminary result for 
fg^^- See reference [11] for a complete set of re- 
sults for decay constants, binding energies and 
wave functions for static- light mesons. For exam- 
ple, for the static-light meson we find 



1.22+; 



(statistical) 



(2) 



fstat 
fstat 

Here we present a sampling of the compari- 
son of the properties of static- light, NRQCD-light 
and NRQCD-NRQCD mesons. Figure 1 displays 
an effective- "mass" plot for the local-smeared 
channel computed using an optimal ("Z") source 
designed to project out the ground state of the 
static-light pseudoscalar meson. Figures 2 and 3 
are similar plots for the NRQCD-light meson and 
the NRQCD-NRQCD meson (bottomonium), re- 
spectively. (The values of the effective masses 
are renormalized differently for different actions 
and so should not be compared.) In all three 
cases, MOST is successful in producing effective- 



mass plots which plateau early, are flat indicating 
ground-state saturation and have small statisti- 
cal errors. Figures 4, 5, and 6 show 15 and 25 
Bethe-Salpeter (Coulomb gauge) wave functions 
(normalized to unity at f = 0) for static-light, 
NRQCD-Hght, and NRQCD-NRQCD mesons, re- 
spectively. A simple NRQCD action (tadpole- 
improved, but with only the kinetic energy term 
with bare lattice mass Ma = 1.7) is used here in 
this illustrative example. Statistical errors seem 
to be as small as in the static case, and will permit 
a precise comparison of the static /s results with 
those forthcoming for an NRQCD heavy quark in- 
cluding 1/M terms in the action and current [12]. 
We plan to apply these techniques to sim- 
ulations using the Mandula-Ogilvie [13] lattice 
implementation of the Isgur-Wise flxed-velocity 
action, where there have been problems in ex- 
tracting a consistent signal using simple smearing 
methods. 
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Static-light, LS, Z-source 



Static-light wave function (measured at t=4) from a Z-source 
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Kappa = 0.154 
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Figure 1. Ground state static- light meson efTec- 
tive mass. 



Figure 4. IS and 2S wave functions normalized 
to unity at the origin for a static-light meson. 



NRQCD-lighl. LS. Z-.source 



NRQCD-light wave function (mea.sured at t=4) from a Z-.source 



20"3 • 30 , bela=6.0. 32 cfg. 
NRQCD: K.E. only, simple update 
n=3. Ma=1.7. u_0=0.88 
Kappa = 0.154 
Z[lS;n#=3;t_0=2.t=3] 
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Figure 2. Ground state NRQCD-Hght meson ef- 
fective mass. 



Figure 5. IS and 2S wave functions normalized 
to unity at the origin for a NRQCD-light meson. 



NRQCD-NRQCD. LS. Z-source 

20'>3 * 30. beta=6.0. 32 cfg. 
NRQCD: K.E. only, simple update 
n=3. Ma=1.7. u_0=0.88 
Z[lS;n#=5;t_0=2,t=3J 





NRQCD-NRQCD 


wave function (measured at t=4) from a Z-.source 
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Figure 3. Ground state NRQCD-NRQCD meson 
effective mass. 



Figure 6. IS and 2S wave functions normalized to 
unity at the origin for a NRQCD-NRQCD meson. 



